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Abstract 

We consider the evolution of a flat, isotropic and homogeneous Friedmann-Robertson- Walker 
Universe, filled with a causal bulk viscous cosmological fluid, that can be characterized by an 
ultra-relativistic equation of state and bulk viscosity coefficient obtained from recent lattice QCD 



calculations. The basic equation for the Hubble parameter is derived under the assumption that the 
total energy in the Universe is conserved. By assuming a power law dependence of bulk viscosity 

> , 

'sj" , coefficient, temperature and relaxation time on energy density, an approximate solution of the 

T— I ! 

OO ■ field equations has been obtained, in which we utilized equations of state from recent lattice QCD 

CN ; 

simulations QCD and heavy-ion collisions to derive an evolution equation. In this treatment for the 

o . 

, viscous cosmology, we found no evidence for singularity. For example, both Hubble parameter and 

scale factor are finite at t = 0, t is the comoving time. Furthermore, their time evolution essentially 



^ ■ differs from the one associated with non-viscous and ideal gas. Also thermodynamic quantities, 

like temperature, energy density and bulk pressure remain finite as well. In order to prove that 
the free parameter in our model does influence the final results, qualitatively, we checked out other 
particular solutions. 

PACS numbers: 04.50.Kd, 04.70.Bw, 97.10. Gz 
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I. INTRODUCTION 



The dissipative effects, including both bulk and shear viscosity, are supposed to play a 
very important role in the early evolution of the Universe. The first attempts at creating 
a theory of relativistic fluids were those of Eckart [ll and Landau and Lifshitz j2|. These 
theories are now known to be pathological in several respects. Regardless of the choice of 
equation of state, all equilibrium states in these theories are unstable and in addition signals 
may be propagated through the fluid at velocities exceeding the speed of light. These 
problems arise due to the first order nature of the theory, that is, it considers only first- 
order deviations from the equilibrium leading to parabolic differential equations, hence to 
infinite speeds of propagation for heat flow and viscosity, in contradiction with the principle 
of causality. Conventional theory is thus applicable only to phenomena which are quasi- 
stationary, i.e. slowly varying on space and time scales characterized by mean free path and 
mean collision time. 

A relativistic second-order theory was found by Israel jsl and developed by Israel and 
Stewart j4|, Hiscock and Lindblom jsl and Hiscock and Salmonson jol into what is called 
"transient" or "extended" irreversible thermodynamics. In this model deviations from equi- 
librium (bulk stress, heat flow and shear stress) are treated as independent dynamical vari- 
ables, leading to a total of 14 dynamical fluid variables to be determined. For general reviews 
on causal thermodynamics and its role in relativity see 

Causal bulk viscous thermodynamics has been extensively used for describing the dy- 
namics and evolution of the early Universe or in an astrophysical context. But due to the 
complicated character of the evolution equations, very few exact cosmological solutions of 
the gravitational field equations are known in the framework of the full causal theory. For 
a homogeneous Universe filled with a full causal viscous fluid source obeying the relation 
^ ~ p^^^, with p the energy density of the cosmological fluid, exact general solutions of the 
field equations have been obtained in js iSj. It has also been proposed that causal bulk 
viscous thermodynamics can model on a phenomenological level matter creation in the early 
Universe [st. Exact causal viscous cosmologies with ^~p'^,S7^1/2 have been considered in 
Ref. 3. 

Because of technical reasons, most investigations of dissipative causal cosmologies have 
assumed Friedmann- Robert son- Walker (FRW) symmetry (i.e. homogeneity and isotropy) 
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or small perturbations around it 12| . The Einstein field equations for homogeneous models 
with dissipative fluids can be decoupled and therefore are reduced to an autonomous system 
of first order ordinary differential equations, which can be analyzed qualitatively jl^ . 

The role of a transient bulk viscosity in a FRW space-time with decaying vacuum has been 
discussed in [U]. Models with causal bulk viscous cosmological fluid have been considered 
recently [isl. They obtained both power-law and inflationary solutions, with the gravita- 
tional constant an increasing function of time. The dynamics of a viscous cosmological fluids 
in the generalized Randall-Sundrum model for an isotropic brane were considered in 16 1. 
The renormalization group method was applied to the study of ho mog eneous and flat FRW 
Universes, filled with a causal bulk viscous cosmological fluid, in 17|. A generalization of 
the Chaplygin gas model, by assuming the presence of a bulk viscous type dissipative term 
in the effective thermodynamic pressure of the gas, was investigated recently in |l8|. 

Recent RHIC results give a stron g in dication that in the heavy-ion collisions experiments, 
a hot dense matter can be formed 191]. Such an experimental evidence might agree with 
the "new state of matter" as predicted in the Lattice QCD simulations [2^. However, the 
experimentally observed elliptic flow in peripheral heavy-ion collisions seems to indicate that 
a thermalized collective QCD matter has been produced. In a addition to that, the success 



of ideal fluid dynamics in explaining several 
spectra of identified particles, elliptic flow 



experimental data e.g. transverse momentum 



2l|, together with the string theory motivated 



that the shear viscosity rj to the entropy s would have the lower limit ^ l/An [22j] leading 
to a paradigm that in heavy- ion collisions, that a nearly perfect fluid likely be created and 
the quarks and gluons likely go through relatively rapid equilibrium characterized with a 
thermalization time less than 1 fm/c 23j. 



According to recent lattice QCD simulations 



the bulk viscosity ^ is not negligible 



near the QCD critical temperature T^. It has been shown that the bulk and shear viscosity 
at high temperature T and weak coupling as, ^ ~ a^T^ / \na~^ and t] ~ T^/(a^lna7^) [26|. 
Such a behavior obviously reflects the fact that near Tc QCD is far from being conformal. 
But at high T, QCD approaches conformal invariance, which can be indicated by low trace 



anomaly (e — 3p)/T'^ 27|, where e and p are energy and pressure density, respectively. In 
the quenched lattice QCD, the ratio (/s seems to diverge near 28|. 

To avoid the mathematical difficulties accompanied with the Abel second type non- 
homogeneous and non-linear differential equations 29|], one used to model the cosmological 
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fluid as an ideal (non- viscous) fluid. No doubt that the viscous treatment of the cosmological 



background should have many essential consequences [30|. The thermodynamical ones, for 
instance, can profoundly modify the dynamics and configurations of the whole cosmological 
background [31]. The reason is obvious. The bulk viscosity is to be expressed as a func- 
tion of the Universe energy density p [32j. Much progress has been achieved in relativistic 
thermodynamics of dissipative fluids. The pioneering theories of Eckart [l| and Landau and 
Lifshitz suffer from lake of causality constrains. The currently used theory is the Israel 
and Stewart theory [3|, , in which the causality is conserved and theory itself seems to be 
stable {5, 7|. 

In this article, we aim to investigate the effects that bulk viscosity has on the Early 
Universe. We consider a background corresponding to a FRW model filled with ultra- 
relativistic viscous matter, whose bulk viscosity and equation of state have been deduced 
from recent heavy-ion collisions experiments and lattice QCD simulations. 

The present paper is organized as follows. The basic equations of the model are written 
down in Section [Tll In Section IIIII we present an approximate solution of the evolution 
equation. Section IIVI is devoted to one particular solution, in which we assume that H = 
const. The results and conclusions are given in Sections |VT] and IVIIt respectively. 



II. EVOLUTION EQUATIONS 

We assume that geometry of the early Universe is filled with a bulk viscous cosmological 
fluid, which can be described by a spatially flat FRW type metric given by 

ds^ = de - a" it) [dr^ + [dO^ + sin^ OdcP^)] . (1) 

The Einstein gravitational field equations are: 

Rik — -gikR = SnGTik- (2) 

In rest of this article, we take into consideration natural units, i.e., c = 1, for instance. 

The energy-momentum tensor of the bulk viscous cosmological fluid filling the very early 
Universe is given by 

Tt = {p + p + n)u,u'-ip + U)St (3) 

where i, k takes 0, 1, 2, 3, p is the mass density, p the thermodynamic pressure, 11 the bulk 
viscous pressure and Ui the four velocity satisfying the condition MjM* = 1. The particle and 
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entropy fluxes are defined according to A^* = nw* and 5** = sN^ — (Tn^/2^T) -u*, where n 
is the number density, s the specific entropy, T > the temperature, ^ the bulk viscosity 
coefficient, and r > the relaxation coefficient for transient bulk viscous effect (i.e. the 
relaxation time), respectively. 

The evolution of the cosmological fluid is subject to the dynamical laws of particle number 
conservation N^.^ = and Gibbs' equation Tdp = d (p/n) + pd (l/n). In the following we 
shall also suppose that the energy-momentum tensor of the cosmological fluid is conserved, 
that is = 0. 

The bulk viscous effects can be generally described by means of an effective pressure 11, 
formally included in the effective thermodynamic pressure Pe// = p + H Tj]. Then in the 
comoving frame the energy momentum tensor has the components Tq = p, = T2 = T| = 
—Peff- For the line element given by Eq. ([1]), the Einstein fleld equations read 

a) =T'^'- « 

d 4:TT ^ , , , , 

- = -yG (3pe// + p), (5) 

where one dot denotes derivative with respect to the time t, G is the gravitational constant 
and a is the scale factor. 

Assuming that the total matter content of the Universe is conserved, T-.j = 0, the energy 
density of the cosmic matter fulfllls the conservation law: 

p + 3Hipeff + p) = 0, (6) 

where we introduced the Hubble parameter H = d/a. In presence of bulk viscous stress 
n, the effective thermodynamic pressure term becomes Pe// = p + H. Then Eq. (E]) can be 
written as 

p + 3H{p + p) = -3UH. (7) 
For the evolution of the bulk viscous pressure we adopt the causal evolution equation 



obtained in the simplest way (linear in 11) to satis 
production to be non-negative, S** = II^/^T > 



the if-theorem (i.e., for the entropy 
). According to the causal relativistic 



Israel- Stewart theory, the evolution equation of the bulk viscous pressure reads \l\ 

Ttl + U = -3^H - irn I 3if + - - I - ^ I . (8) 



In order to have a closed system from equations (jlj), ([7]) and ([8]) we have to add the equations 
of state for p and T. 

As shown in Appendix A, the equation of state, the temperature and the bulk viscosit y of 
the quark-gluon plasma (QGP), can be determined approximately at high temperatures 25| 
from recent lattice QCD calculations js^, as 

P = ojp, T = Pp\ ^ = ap + —T^, (9) 

with w = (7 - 1), 7 ~ 1.183, r ~ 0.213, f3 ~ 0.718, 

1 97^-247 + 16 



a 



(10) 



9uo 7 — 1 

and Uq ~ 0.5 — 1.5 GeV. In the following we assume that ap » d/ujoT^, and therefore we 
take ^ ~ ap. In order to close the system of the cosmological equations, we have also to 
give the expression of the relaxation time r, for which we adopt the expression j?!, 

r = ^p~^ ~ a. (11) 

Eqs. ([9]) are standard in the study of the viscous cosmological models, whereas the equa- 
tion for r is a simple procedure to ensure that the speed of viscous pulses does not exceed 
the speed of light. Eq. ffTTl) implies that the relaxation time in our treatment is constant 
but strongly depends on EoS. These equations are without sufficient thermodynamical mo- 
tivation, but in the absence of better alternatives, we shall follow the practice of adopting 
them in the hope that they will at least provide some indication of the range of bulk viscous 
effects. The temperature law is the simplest law guaranteeing positive heat capacity. 

With the use of Eqs. ([HD, dS]) and (ITTi) . respectively, we obtain the following equation 
describing the cosmological evolution of the Hubble function H 

H+l[l + (l- r)^]HH +-H-(1 + r)H-'H^ + ?(7 - 2)i/' + ^-H^ = 0. (12) 
2 a Ala 

III. AN APPROXIMATE SOLUTION 

We introduce the transformation u = H, so that Eq. fll2p is transformed into a first order 
ordinary differential equation, 

du , , 10/3, , A 91o37o 



""dH'^^^ + (,2 + ~ + " J ^ + 4 (^^ + 2 a = °- ^^^^ 
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We can rewrite Eq. (fT3|) in the form 



where 



n^=F^{H)n + Fo{H), (14) 



1 + r 

Qs = u E = u exp ( — / dH 



H 



Fo{H) = -i^-{^-2)H' + --^H'jE\ 
By introducing a new independent variable z = J Fi{H) dH, we obtain 



n^-n = g{z), (15) 



with g{z) is defined parametrically as, 



9i^) = (16) 



As shown in Appendix B, g{z) can be approximated as a simple function of z 

g{z) ^Cz, (17) 

where C is a constant. We proceed with this approximation to get solvable differential equa- 
tions. Keeping the parametric solution of g{z), Eq. (ISUj) . results in much more complicated 
differential equations. This would be the subject of a future work. 
From the definitions of Vt and z we have 

VL = H^+'H, (18) 

V 2(l-r) ar J ' ^ ^ 

Analogous to the solution of reduced Abel type canonical equation, 

- y = ax (20) 

(see Appendix C) we obtain the relation Q = z/V. Therefore, from Eqs. f|T8l) and f|T9|) we 
obtain the following first order differential equation for Hubble parameter H, 

VH = Z*±(^lMff^ + Lh (21) 

2(1 -r) ar ^ ' 



with the solution 

»m = ,,p(_af;y)_^ (22) 



where 



^.zfcl^lM. (23) 

2 1 — r ar 



and P is taken as a free parameter. We can assign any real value to V. For the results 
presented in this work, we used a negative value. This negative sign is necessarily to overcome 
the sign from the integral limits. The geometric and thermodynamic quantities of the 
Universe read 

aW=a„f ""'-f,^' Y'\ (24) 

^\exp{-Bt/V)-Aj ' ^ ^ 

^W = '^''^=K^ [expW)--4|0 ^- 

^ ^ ' p V [exp(-5t/p) - y ' ^ ^ 

g(t) = - 1 = exp(-fit/P) - 1. (28) 

Oq is an arbitrary constant of the integration. The sign of q indicates whether the Universe 
decelerates (positive) or accelerates (negative), q can also be given as a function of the ther- 
modynamic, gravitational and cosmological quantities g(t) = [p{t)+3p(t) + 3Il{t)]/2p(t) |33 |. 



IV. DE SITTER UNIVERSE 

Besides the approximation in g{z), previous solution apparently depends on the free 
parameter V. In this section, we suggest a particular solution to overcome V. Eq. fll2p can 
easily be obtained by assuming that H doesn't depend one t, i.e, de Sitter Universe. With 
a simple calculation, we get an estimation for H 

H=^^. (29) 
92-7 ^ ^ 
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The geometric and thermodynamic parameters of the Universe are given by 



a{t) = ooexp 



9(2-7) 



Pit) 



T(t) = 3"/3 



[9(2-7) 
4a"^7 



n(t) = -37 

g(t) = -1. 



[9(2-7) 
4a~"^7 
9(2-7) 



2r 



-\ 2 



(30) 

(31) 

(32) 

(33) 
(34) 



Although we have assumed here that the cosmic background is filled with viscous matter, 
the assumption that H = const results in an exponential scale parameter, Eq. f pOj) . This 
behavior characterizes the de Sitter space, when A = A; = 0. p and T are finite at small t as 
given in Fig. |2j 



V. PARTICULAR SOLUTION 

Another particular solution for Eq. (fT2|) can be obtained, when assuming that the depen- 
dence of u on H can be given by the polynomial in Eq. ( 12T1) 



u = biH^ + b2H, (35) 

where 61 and 62 are constants. Some simple calculations show that this form is a solution of 
the initial equation, Eq. (fT3l) . if 

u 31 + 7 

h = (36) 

62 = — . (37) 

ra 

62 is identical to B in Eq. fl23|) . r and 7 have to satisfy the compatibility relation 

2-7 
r = . 

2 + 72 

Integrating Eq. fl35|) results in 



(38) 



1 - 61 exp(-02t) 
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where minus sign in the exponential function refers to flipping the integral limits. This was 
not necessary while deriving the expressions given in Section IlIII The free parameter P 
compensates it. The geometric and thermodynamic quantities of the Universe read 




2r 



(40) 



(42) 



(41) 



^^^^ ^ bl [2exp(M)-37] 

[exp{b2t) - bif 
q{t) = exp{b2t) - 1. 



(43) 



(44) 



Obviously , we notice that the scale parameter in Eq. fHOj) looks like Eq. 1^^, which strongly 
depends on the free parameter V. The other geometric and thermodynamic quantities find 
similarities in Eq. ( l25l) - ( 128|) . respectively. Deceleration parameter q seems to be positive 
everywhere. 

VI. RESULTS 

In present work, we have considered the evolution of a full causal bulk viscous flat, 
isotropic and homogeneous Universe with bulk viscosity parameters and equation of state 
taken from recent lattice QCD data and heavy-ion collisions. Three classes of solutions of 
the evolution equation have been obtained. 

In Fig. [H H{t) and a{t) are depicted in dependence on the comoving time t. We compare 
H{t), given by Eq. f l22|) . and a{t), given by Eq. fl2^ . with the counterpart parameters 
obtained in the case when the background matter is assumed to be an ideal and non- viscous 
fluid, described by the equations of state of the non-interacting ideal gas. 



In the left panel of Fig. [H H{t) = a/a has an exponential decay, whereas in the non- 
viscous case, H{t) is decreasing according to Eq. (H5|) . The latter is much slower than 
the former, reflecting the nature of the exponential and linear dependencies. The other 



H{t) 



1 



(45) 
(46) 



2t 



a 



(t) = Vt. 
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Fig. 1: Left panel: H{t) vs. t for viscous (solid line) and non-viscous fluid (dashed line). Contrary 
to non-viscous fluid, H{t) is our solution shows no singularity. The straight line gives the results 
from the particular solution, Eq. (j29p . Double dotted curve represent the particular solution, 
Eq. (|39|) . Right panel: a{t) vs. t. The approximate solution gives finite a{t) at i = 0. (solid line). 
Dotted straight line represents the results from Eq. (j30p . The particular solution is given by the 
double dotted hue, Eq. (|iO]) . 



IS nuid 
non-viscous fluid 
particular 





Fig. 2: Energy density p{t) as a function of t (left panel). The dependence of T on t is given in 
the right panel. The curves are as in Fig. [TJ In both case, viscous fluid gives no singularity at 
vanishing t. Straight curve are from Eq. (j3ip and (j32p . respectively. 



difference between the two cases is obvious at small t. We notice a divergence, or singularity, 
associated with the ideal non- viscous fluid, Eq. (I45l) . The viscous fluid results in finite H 
even at vanishing t, as can be seen from Eq. fl2^ . 

The scale factor a{t) also shows differences in both cases. a{t) in a Universe with an 
ideal and non- viscous background matter depends on t according to Eq. (146|) . which simply 
implies that a{t) is directly proportional to t, and a{t) vanishes at t = 0, which shows the 
existence of a singularity of H. Assuming that the background matter is described by a 
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t t 

Fig. 3: Left panel: bulk pressure 11 depicted in dependence on t. At small t, U jumps from 
negative to positive values. At larger t, 11 vanishes as the case in the non-viscous fluid (dashed 
line). The bottom line gives 11 according to first particular solution, Eq. ()33p . Right panel: 
deceleration parameter q is depicted with t. Straight and dotted lines represent non-viscous fluid 
and first particular solution, Eq. (|34|l . respectively. Top curve represents second particular solution, 
Eq. dSD. 

viscous fluid results in different a(t)-behaviors with increasing t. At t = 0, a{t) remains 
finite. Correspondingly, H{t) remains also finite. In general, the dependence on t is much 
more complicated than in Eq. (146|) . Here we have an A/V root of an exponential function. 
If exp{—Bt/V) » A, a remains constant. 

Fig. |2] illustrates the dependence of the two thermodynamical quantities, p and T, on 
the comoving time. The non-viscous Universe shows a singular behavior in p at vanishing 
t, as shown in the left panel of Fig. [2j This is not obvious in the case where we have taken 
into consideration a finite viscosity coefficient, i.e., p is finite at t = 0. In both cases, p is 
decreasing with increasing t, reflecting that the Early Universe was likely expanding. Also 
the life time of the thermal viscous Universe seems to be shorter than for the non-viscous 
Universe. Almost the same behavior is observed in the right panel of Fig. [21 The temperature 
T seems to be finite at vanishing t in the viscous Universe. The T-singularity is only present, 
if we assume that the background matter is non-viscous ideal gas. 

In left panel of Fig. [3], we show the dependence of the bulk viscous pressure 11 on t. 11 
takes negative values at very small t. Then it switches to positive values at some values 
of t. After reaching the maximum value, 11 decays exponentially with increasing t. At 
larger t, U entirely vanishes. The deceleration parameter q, given by Eq. (l28l) . is depicted 
in the right panel of Fig. [3l and it is compared with q for a non-viscous fluid, q = —3. 
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The approximate solution, given by Eq. (l28i) . results in negative q at small t, referring to 
expansion era. q from the particular solution, Eq. flM|) is negative everywhere. 

For the particular solution, only the scale factor depends on t, Eq. fl5U]) . The results are 
given in the right panel of Fig. [T] All cosmological and thermodynamical quantities given 
by Eq. (l29l) and Eqs. (!3T!) - (!34l) are constant in time. 

VII. CONCLUSIONS 

It is obvious that the bulk viscosity plays an important role in the evolution of the Early 
Universe. Despite of the simplicity of our model, it shows that a better understanding of the 
dynamics of our Universe is only accessible, if we use reliable equation of state to characterize 
the matter filling out the cosmic background. 

We conclude that the causal bulk viscous Universe described by the approximate solution 
starts its evolution from an initial non-singular state with a non-zero initial value of Hubble 
parameter H{t) and scale factor a{t), where t is the comoving time. In this treatment, t 
is given in GeV~^. Also the thermodynamical quantities, energy density p for instance, 
are finite at vanishing t. Even the temperature T itself shows no singularity at t = 0. The 
Hubble parameter H decreases monotonically with T similar to p. The bulk viscous pressure 
n likely satisfies the condition that H < at very small t indicating to inflationary era. Then 
H switches to positive value. It reaches a maximum value and then decays and vanishes, 
exponentially, at large t. The deceleration parameter q shows an expanding behavior in the 
case of non-viscous ideal gas and first particular solution. For second particular solution, q 
starts from zero and increases, exponentially. According to this solution, the Universe was 
decelerating. The approximate solution shows an interesting behavior in q{t), Eq. ( |28i) . At 
small t, the values of q are negative, i.e. the Universe was accelerating (expansion). At larger 
t, a non- inflationary behavior sets on, g > 0, i.e., the Universe switched to a decelerating 
evolution. 

In this treatment, we assumed that the Universe is flat, k = 0, and the background 
geometry is filled out with QCD matter (QGP) with a finite viscosity coefficient. The 
resulting Universe is obviously characterized by a shortly increasing and afterward constant 
scale factor and a fast vanishing Hubble parameter. At t = 0, both a{t) and H{t) remain 
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finite, i.e., there is no singularity. The vahdity of our treatment depends on the validity of 
the equations of states, Eq. [9l which we have deduced from the lattice QCD simulations at 
temperatures larger than Tc ~ 0.19 GeV. Below T^, as the Universe cooled down, not only 



the degrees of freedom suddenly increase 



36l ] but also the equations of state turn to be the 



ones characterizing the hadronic matter. Such a phase transition - from QGP to hadronic 
matter - would characterize one end of the validity of our treatment. The other limitation 
is the very high temperatures (energies), at which the strong coupling as entirely vanishes. 
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Appendix A: Viscosity coefficient ^{T) from LQCD 



Following the discussion presented in [3J|, the bulk viscosity of QGP can be calculated 
from the lattice QCD by Eq. (13) in that paper. We assume that the decay factors for pions 
and kaons are vanishing above the critical temperature of the phase transition QGP-hadrons. 
The quark- ant iquark condensates can be neglected at temperatures higher than the critical 



one f35|. Therefore, Eq. (22) of Ref. [3J] would be reduced to 



9cc;oe = Ts (^^-3^ -4(p-3p) + 16|e,| (47) 

where p is the energy density and = dp/ dp is the square of the speed of sound. The 
parameter Uo is a scale depending on the temperature T, and defines the validity of the 
underlying perturbation theory. In this relation, the viscosity is assumed to have a thermal 
part which can be determined through lattice calculations, and a vacuum contributing part, 
which can be fixed using quark and gluon condensates. The vacuum part would take the 
value 

16|e^|(l + I ■ 1.6) ~ (560 MeV)^ ~ {3TcY (48) 
8 

Our algorithm is the following. Using lattice QCD results on trace anomaly, (e — 3p)/T^, 
and other thermodynamical quantities, we can determine the bulk viscosity. To make use of 
the lattice QCD results, it is useful to make a suitable fit to the data at high temperatures. 
Then we obtain the following equations of state 

p = up, T = c1 = uj 

where u = 0.319, (3 = 0.718 ± 0.054 and r = 0.23 ± 0.196. Using the equations of state, 
Eq. (gg]) in Eq. (gTD, we obtain 

1 -24.^16 (49) 

9u)o 7-1 tUo 

Appendix B: Estimations of g{z) 

For analytical purposes, the function g{z), which is defined in z parameter as g{z) = Fq/Fi 
in Eq. (fT6|) . can be numerically estimated depending on the parameter z by using the 
following procedure. First, we plot it parametrically depending on the parameter H, Fig. 
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Fig. 4: The parametric dependence of g{z) on z is given by open symbols. Eq. ()50p is depicted as 
solid line. The dotted line represents the linear fit, Eq. ()5ip . 



dlj. Then we fit the resulting curve to various functions. Based on least-square fit, best 
choice would be a mixture of polynomial and exponential functions, 

exp{dz) + e 



g{z) = a + bz + c- 



(50) 



[expidz) + ff' 

where the coefficients read a = -2.078 ± 0.117, b = 0.091 ± 0.007 and c = 17.332 ± 1.553, 
d = 0.189 ±0.003, e = -0.814 ±0.162 and / = 2.849 ±0.02. At small values of z, it is clear 
that the dependence is linear, 

g{z) = c + Cz. (51) 

Obviously, the intersect c is much smaller than the slope C. The sign of g{z) can be flipped 
regarding to the sign of its independent variable z. Accordingly, we get 



g{z) ^ Cz. 



(52) 



To prove this dependence, algebraically, we try to estimate g{z) directly from the division 
of Fq by -Fi, which can be approximated by including their first terms only, i.e. 

3(7 - 2) 



!>m " 2|i + (1 _ ,)^J 
Then, we approximate z{H) to the form. 



H 



l-r 



2(1 -r) 



(53) 



(54) 
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Finally, we now able to derive an approximate estimation for g{z). According to Eq. (153|) 
and fl54|) . we get 

a(z) - . (55) 

Amazingly, this expression looks the same as the one we obtained from the numerical ap- 
proximation with 

C=(i^lM. (56) 
|l + (l-r)7f 

Appendix C: Solution of Abel equation yy — y = ax 

To solve Eq. fl2U]) we divide the whole equation by and introduce a new variable 
V = 1/y. Then Eq. fl20|) reads 

dv 

— + f ^ + axv^ = 0. (57) 
ax 

We then introduce the function v = w/x. 

dW 2 n 

X—— = w — w — aw , (58) 

ax 

Previous differential equation can be solved by separation of variables 

dw 



w — w'^ — aw^ 



InC'^x, (59) 



where C is an arbitrary constant of integration. To calculate the integral, we write the 
function to be integrated as 

1 1 aw 1 



w — w'^ — aw^ w aw'^ + w ~ 1 aw'^ + w — 1 
Let us assume that A = 1 + 4a > (this implies that a > 0). 



(60) 



/ 



dw 1 , 2aw — \/A + l 1, . 2 i\ 1 /^.N 

In ^ in i^aw +w — l)+mw. [bl) 



w-w^-aw^ 2VA 2aw + VA + l 2 
Therefore the general solution of Eq. (!58l) can be written as 

1/2VA 



W 



Vaw^ + W-1 \2aw - VA + 1 ' 



leading to 

ti' Vaw2 + _ 1 \ 2aw - V A + 1 / 
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